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ABSTRACT 

Both radiative and mechanical feedback from Active Galactic Nuclei have been found to be impor- 
tant for the evolution of elliptical galaxies. We compute how a shock may be driven from a central 
black hole into the gaseous envelope of an elliptical galaxy by such feedback (in the form of nuclear 
winds) using high resolution 1-D hydrodynamical simulations. We calculate the synchrotron emission 
from the electron cosmic rays accelerated by the shocks (not the jets) in the central part of ellip- 
tical galaxies, and we also study the synchrotron spectrum's evolution using the standard diffusive 
shock acceleration mechanism, which is routinely applied to supernova remnants. We find quanti- 
tative consistency between the synchrotron radio emission produced via this mechanism with extant 
observations of elliptical galaxies which are undergoing outbursts. Additionally, we also find that 
synchrotron optical and X-ray emission can co-exist inside elliptical galaxies during a specific evolu- 
tionary phase subsequent to central outbursts. In fact, our calculations predict a peak synchrotron 
luminosity of ~ 1.3 x 10^ Lq at the frequency 5 GHz (radio band), of ^ 1.1 x 10^ Lq at 4.3 x 10^'' 
Hz (corresponding to the absolute magnitude -10.4 in R band), and of 1.5 x 10^ Lq at 2.4 x 10^^ 
Hz (soft X-ray, 0.5 — 2.0 keV band). 

Subject headings: galaxies : ISM — radio continuum : galaxies — ISM : cosmic rays — radiation 
mechanisms : nonthermal — ISM : jets and outflows 



1. INTRODUCTION 

Is there direct evidence for feedback from Active Galac- 
tic Nuclei (hereafter AGN) available from the non- 
thermal emission of the central regions of massive galax- 
ies? It has been observed that 10-20 percent of AGN 
are radio-loud and that the host galaxies of radio- 
loud AGNs are usu ally massive elliptical galaxies (e.g., 
iHooper et all Il995t |Best 2007). The radio emission 
in radio galaxies and quasars comes from a variety 
of morphological structures, with different characteris- 
tic sizes. For example, there are compact radio cores 
(usually with flat power law spectra) coincident with 
the nucleus of the associated optical objects, as well 
as large scale radio lobes. Other systems show radio 
lobes which may extend to several hundred kiloparsecs 
to a few megaparsecs and usually show steep power 
law r adio spectra (e.g., IKellermann et al.l 11994 IKrolikI 
19991: iK cmbhaavi fc Narlikari 119991: IGilbert et al. 1 120041 : 
Mullin et al.i i2006j) . The typical core radio luminosity 
for FR-II sources at 5GHz is ~ 10^*^ erg s'^ Hz'^, 
while the total luminosity at 178 MH z is ~ 10^'^ erg 
s^^ Hz~^ (e.g., iChiaberge et al.l [20001 ) . with the spec- 
tral index (absolute value) typically smaller than 0.5 at 
the very center and ~ 0.7 — 0.8 at a larger radius. In 
addition to the radio observations, there are also obser- 
vations in other bands for radio galaxies. For example, 
in the R-band 3CR Snapshot Survey of 252 radio galax- 
ies , the wind-like stru ctures seen in the contour figures 
of lMartel et all ([l999l ) may be interpreted as produced 
by nuclear winds. Kiloparsec-scale winds are also ob- 
served in high redshift quasars (e.g.,'de Kool et al. 2001'; 
Nesvadba 2009 : Ncsyadba et al. 2008; Moc et al. 2009; 
Cattaneo et al I l2009l and references therein). Optical 



synchrotron jets ar e also confirmed in many sources (e.g., 
IMartel et al.f[T999l) and can extend to distances far out- 
side of the body of the galaxies. Furthermore, ultravio- 
let and X-ray observations have also found hi gh-velocit; 
outflows in some quasars (e.g., Nesvadba et al.l 1200 



Hamann et all I2008t iNesyadbal 120091: lAlexander et all 
20091: IChartas et al 1 120091 1 Finally, although the ev- 
idence is not absolutely conclus ive, there are observa- 
tions (e.g.. IChiaberge et al.ll200dl ) indicating optical syn- 
chrotron from the cores of some active elliptical galax- 
ies. The obvious question is "what processes accelerate 
cosmic rays (hereafter CR) electrons up to the energies 
required to produce the observed synchrotron emission?" 

Synchrotron emission is believed to be the mecha- 
nism responsible for the radio emis sion in AGNs , based 
on their strong polarization (e.g., IKrolikI ri999( ). The 
relativistic electrons, which are required for the syn- 
chrotron radio emission, can be accelerated by several 
sources. Jets are observed in most of the radio-loud 
AGNs and most investigators believe that th ose jets pro- 
duce this radio synchrotron emission (e.g., fReallgnl; 
iBlandford fcRe es 1974; Sironi & Socrates 2009; also 
see the review [Mirabel fc Rodrfgueail999l and references 
therein). However, the energy carried by the jets is only 
dissipated in small areas near the end point or within in- 
ternal reflected shocks. As the jets are driven outwards 
with relativistic velocity, the centers of the radio con- 
tours from the jets do not overlap with the centers of the 
galax ies (see, e.g., rKellcrmann et al. 1994; MuUin ^t al.l 
I2006D . So, ahhough jets are very likely responsible for 
the displaced radio lobes, there must be other mecha- 
nisms responsible for the observed radio core emission 
at the centers of FR-II galaxies. Furthermore, in or- 
der to explain synchrotron emission in galaxies with- 
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out jets, other sources able to accelerate electrons are 
needed. In fact, jets are observed in ~ 80% of FR-I ra- 
dio galaxies, in ~ 40-70% radio quasars ^ and even fewer 
in highly l uminous FR-II galaxies (e.g.. iBridle &: PerlevI 
[1984; Ke mbhaavi &: Narlikad[T999h . 

This paper addresses the inevitable generation of 
radio synchrotron sources in elliptical galaxies conse- 
quent to central outbursts. In fact, it is now widely 
accepted that supermassive black holes (SMBHs) within 
a mass range of 10^ ^ 10^'^ Mq reside at the center of 
bulges and elliptical galaxies, and that feedback from 
these SMHBs can profoundly a ffect the formation and 
evolution of these galaxies f e.g.. iKormendv fc Richstone 
1995; Silk & Rees 1998; F abianT 
Tool; King 2003; Granato ct aLl 



9991 IBurkert fc Silk 
20041: 'Springel et al 



2005; Sazonov et al. 2005; HopkinseLaD'2006; Hu et al 
2006f ). This is strongly supported by the observed rela- 



properties of their host galaxies (e.g.. iMagorrian et al. 


1998 


; IGebhardt et a 


.1 120001: iFerrarese & Merritt 


200C 


; iTremaine et al.l 


20021: McLure & Dunlod 120021: 


Gral 
2008 


lam & Driveii 12007 


iLauer et al. 20071: iLou & Jiang 


). In particular, in a series of works based on high 



resolution numerical simulations whic h include both 
radiative and m echanical feedback (jCiotti &: OstrikeJ 
[19971 2001, 2007: lCiotti"eraIll2009l hereafter Paper I), 
it has been shown how the mass loss from evolving stars 
can drive significant nuclear activity, characterized by 
strong and recurrent nuclear bursts, even in the absence 
of merging. This model su pports the idea that the 
"cooling flow" phases (e.g., iPeterson fc Fabianl 12006( 1 
and quasar phases are different aspects of the evolution 
of a normal elliptical galaxy. In the calculations, it is 
found that accretion occurs in bursts, during which both 
radiative and mechanical output from the central SMBH 
pushes matter out and drives shocks into the galactic 
gas. The outbursts typically combine three physically 
separate phenomena: central star bursts, mechanical 
feedback from winds emanating from broad-line regions 
surrounding the SMBH, and radiative feedback from 
absorption and scattering of hard X-ray photons. The 
resulting shocks are similar in kind (but more modest 
in degree) to the predicted outflowing blast waves 
produced by AGNs after galaxy merging (e.g ., see 
Di Mattc o et al.l l2005t [S oringcl et al. 2005; Naa b et all 
2006t iSiiacki et al.ll2007t see also Johansson et al.ll2009i r 



In this paper, we focus attention on the observa- 
tional properties of the emitted synchrotron emission 
from the shocks that result from the recurrent bursts. 
In fact, since the computed outflow velocities (up to 
several thousand kilometers per second) and densities 
(0.1 to 10 particles per cubic centimeter) are similar to 
what is obser ved in Galactic SNRs such a s Tycho and 
Kepl er (e.g., iCassam-Chenai' et al.l l2007t iDickel et al.l 
one would expect that the same processes as 
those we see acting in SNRs on the pc scale (specif- 
ically at the shocks' surfaces, Wang 2008) would act 
on the kpc scale in galaxies. One can expect there 
to be efficient acceleration of both ionic and electronic 
cosmic rays via the well-known diffusive shock accel- 
eration mechanism (also known a s first-order Fermi 
shock acceleration, e.g., iBlandfo rd fc Ostriked 119781 : 
IBlandford fc Eichled [19871 : iBerezhkol: Enisonlll999r see 



also iTreumann fc Jaroschekl l2008f) . Then, just as the 
accelerated CR electrons produce synchrotron emission 
through standard interactions with the (possibly ampli- 
fied) magnetic field within normal SNRs, so we should 
expect very similar processes to occur in the ISM of el- 
liptical galaxies after outbursts. To our knowledge, this 
has not yet been studied in quantitative detail. Using 
diffusive shock acceleration theory, we can make a rough 
estimate of the scaling relation for luminosity, because 
the synchrotron luminosity only depends on the CR elec- 
tron density, the magnetic field B and the shock velocity. 
If the density and velocity of the shock are fixed, the lu- 
minosity is determined by the volume occupied by CR 
electrons and the magnetic field. The extent of emitting 
region near the shock is determined by the synchrotron 
losses of the electrons, which are proportional to B^, 
so that the volume occupied by CR electrons goes with 

/ B^ (r is the radius of the shock). Because power per 
volume radiated by the electrons goes like B^, the total 
luminosity scales approximately as r^. As the size of an 
elliptical galaxy is ^ 100 times larger than the SNRs, 
and the densities and shock velocities in the cores of el- 
liptical galaxies are expected to be very similar to those 
in SNRs, we expect the synchrotron luminosity from el- 
liptical galaxy cores to be ~ 10'* — 10^ larger than the 
luminosity from SNRs. Our detailed calculations (see §3) 
confirm these rough scaling estimates. 

In this paper, we calculate synchrotron emission from 
the relativistic electrons accelerated by shocks during the 
evolution of an elliptical galaxy and we also show how 
the synchrotron spectrum evolves with time. The Ciotti 
fc Ostriker model for the evolution of elliptical galaxies 
provides us with a specific scenario for computing the 
shocks, which are formed repeatedly due to the feed- 
back of the central SMBHs. With the propagation of 
the shock into the interstellar medium (ISM), electrons 
will be accelerated at different positions (via the stan- 
dard diffusive shock acceleration mechanism), and high 
energy particles can be found at large radii. This is very 
different from the concept sometimes proposed that CRs 
are accelerated at the center and escape to large radius, 
since such CR electrons would suffer from large adiabatic 
losses and become relatively ineffective radiators. In §2 
we present the method to calculate the accelerated elec- 
tron spectrum and associated synchrotron emission. In 
§3 we show the result from a single shock taken from 
a numerical simulation. The main results are summa- 
rized in §4 and discussions and observational results of 
our model are given in §5. 

2. CALCULATION METHOD 

In this Section, we describe the numerical model 
adopted for the simulation and we give the formulae 
for the diffusive shock acceleration mechanism and syn- 
chrotron emission used in this paper. 

2.1. The Hydrodynamical Model 

The details of hydrodynamical simulations of radiative 
plus mechanical feedback and how the kinetic energy, 
momentum and mass of the wind are transferred to the 
ISM in combined (i.e., radiative plus mechanical feed- 
back) models are described in Paper I, while the explo- 
ration of parameter space in combined models is post- 
poned to iCiotti et al.. (,2009g) (hereafter Paper HI) . We 
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recall that the code is 1-D, and that, as in Paper I, a sim- 
plified version of mechanical feedback is adopted^, which 
is sim ilar to that used by others (e.g., iDi Matteo et al.l 
120051 ) studying AGN feedback. The galaxy model is a 
Jaffe stellar distribution embedded in a dark halo so that 
the to tal density profile decreases as (ICiotti et al.l 
l2009bl and re ferences therein). T he cooling and heat- 
ing functions (jSazonov et al.|[2005f ) include photoioniza- 
tion plus Compton and line heating, while the SMBH 
accretion rate is mediated by a circumnuclear accretion 
disk modelized at the level of subgrid physics. Finally, 
the mass return rate from evloving stellar populations is 
computed by using the detailed prescriptions of stellar 
evolution. A discussion of the limitations of 1-D simula- 
tion is deferred to 21 

In the following discussion, we focus on a specific 
model labeled B3o2 (Paper I; Paper III, in preparation), 
which models the radiation and mechanical feedback 
from AGNs to the host elliptical galaxies. This model 
has an initial stellar mass Mj,=2.87 x 10^^ M© and an 
initial SMBH mass Mbh=2.87 x 10^ Mq. The galaxy 
effective radius is Re=6.9 kpc and the aperture central 
velocity dispersion is cro=260 km s""'^. The first grid point 
in the code corresponds to i?inin =5 pc. The simulation 
starts at 2 Gyr and lasts to 14.5 Gyr. The time evolution 
of the SMBH bolometric accretion luminosity is shown 
in Fig. m 

The com putations show that a cooling instability 
(|Fieldlll965D starts as a cool (T - lO'^ K) and dense shcU 
at r 1 kpc; this shell collapses, feeds gas into the center, 
and the SMBH responds with a major burst. In more de- 
tail, the infalling shell compresses the gas inside the col- 
lapsing volume, so that more and more gas is accreted 
by the SMBH even before the cold shell reaches the cen- 
ter. The luminosity of SMBH increases rapidly with 
time, and when it reaches ~1 % of Eddington luminosity, 
the infalling ma tter is pushed out by the radiation from 
the S MBH fe.g.. lOstriker et al.|[l976l : iMilosavlievic et all 
l2009f) . At the same time, a wind appears from the cir- 
cumnuclear disk and a shock is driven into the ISM. The 
shock hits the still infalling cold shell at a very high Mach 
number, due to the low temperature of the shell. The 
shock continues to propagate into the galaxy; gradually 
its energy is deposited into the ISM, and it dies away with 
time. In the meantime, reverse shock waves carry fresh 
gas to the center, and several sub-bursts are generated. 
Finally, the last burst ends the sequence. The duty cycle 
— the fraction of time that the central SMBH is in the 
outburst mode with Lbh > LEdd/30 — is of the order of 
~ 5 %, consistent with the fraction of elliptical galaxies 
seen as quasars. As commonly found in combined (i.e., 
radiative plus mechanical feedback) models discussed in 
full extent in Paper HI, the temporal structure of each 
major burst is highly organized into a series of sub-bursts 
(as apparent from the temporal structure of each of the 
4 major bursts in Fig. [1]), while purely mechanical feed- 
back models present much simpler bursts, as can be seen 
by comparison with Figs. 2, 3 and 4 in Paper I. In our 
case, each shock will accelerate electrons and protons. 

Now we focus on the first sub-burst in the last ma- 
jor burst around 6.5 Gyr. The SMBH bolometric accre- 

^ In practice, the time derivative in equation (28) of Paper I is 
set to zero. 




t(Gyr) 



Fig. 1. — Time evolution of the SMBH bolometric accretion lu- 
minosity Lbh fof model B3o2. The nearly horizontal solid line is 
the Eddington luminosity LEdd- There are four major bursts dur- 
ing the evolution history of the model, and each major burst is 
composed of several sub-bursts. 

tion luminosity Lbh reaches 1 % of Eddington luminosity 
LEdd at the time 6.48 Gyr, and we define this time as the 
zero time point. With this time origin, Lbh reaches 0.1 
LEdd around 0.234 Myr later, and it peaks at 0.5LEdd 
at At = 0.252 Myr. Fig. [2] shows the radial profile of 
the physical variables at At = 0.204 Myr, just before 
the peak of the burst. We can see clearly the cold shell 
initially at 40 pc falling to the center. Fig. [3] shows the 
appearance of outflow with the solid line {At = 0.235 
Myr) still before the peak of the burst: the wind pushed 
out by radiation from the SMBH reaches 5.1 pc, with 
a velocity of only 100 km/s (arrow A). By the time the 
shock reaches 15 pc (arrow B) the maximum velocity has 
increased up to 3600 km/s, as it is being steadily accel- 
erated from below, and it has hit the infalling cold shell, 
which results in a very high Mach number shock. 

2.2. Acceleration of Particles at the Shocks 

Synchrotron emission from shocks based on the dif- 
fusive acceleration mechanism has been applied to as- 
trophysical systems on different scales. In particular, for 
SNRs there exist different nonlinear models for the struc- 
ture and evolution of cosmic ray modified (strong) shocks 
(e.g., see Kovama et al. 1995; Bcrczhko & Ellison 199^ 
EUison ct al. 2000; Lazcndic ct al. 2004; El lison et all 
,2004 .Berezhko fc Volk. 2006: Ellison et al..,200" 7l: see the 
reviews bv lRevnoldj2008llTreumann fc Jarosch ck 2008). 
With this method, the synchrotron spectra at different 
radii and their evolution with time can be computed and 
then compared with the observations of SNRs, show- 



ing g o od agreement (e.g., .Reynolds 1998; .Lazendic et al.l 
120041 : iCassam-Chenai' et al.l 120071 ) . Some researchers 
have also carried out simulations on the scale of galaxy 
clusters to study the radio haloes and relics based on the 
diffusive shock acceleration mechanism in the limit of test 
particle theory (e.g.. .EnBlin ct al...2007 : Pfrommcr et al.l 
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Fig. 2. — Thermal instability and the falling cold shell: the radial 
profile of the ISM energy per unit volume E, proton number density 
Up, radial velocity v, and temperature T. Solid, dotted and dashed 
lines correspond to At = 0.204, 0.214 and 0.224 Myr after Lbh 
reached 0.01 Lecu (6.48 Gyr, see Fig. [ij. This time is defined to 
be the zero point. Note how the gas density in the central regions 
increases due to compression before the cold shell (40 pc from the 
center at At = 0.204 Myr) reaches the center. 

I2008f l. Here, we apply the same techniques to the case 
of wind-driven shocks in elliptical galaxies. 

Our CR shock acceleration scheme is based on well es- 
tablished physical principles (the first order Fermi mech- 
anism) but utilizes a new, quite efficient numerical im- 
plementation, which allows us to perform the analysis 
during post-processing instead of during the runs, sav- 
ing considerable computational time. We use the test 
particle theory to calculate the electron spectrum at 
the shock, which gives the basic characteristics of the 
non-thermal electro ns (e.R ., iBlandford fc Ostrikeii 119781 : 
Blandford & EichldllMa see also the recent review by 
Reynolds 2008i). Consistent with test particle theory, 
the CRs are assumed to not change the structure of the 
shocks, so that we can use the output of the hydrody- 
namical code to calculate the synchrotron spectrum. 

In order to compute the acceleration of the electrons 
and protons, first we have to identify and follow the 
shocks as they move through the ISM. Each shock is 
driven by the combined effect of radiation and nuclear 
winds following a burst, and it propagates into the ambi- 
ent gas of the elliptical galaxy. For each hydrodynamical 
output file, we calculate the dimensionless compression 
factor 

S ^ /^div(v) (1) 

\v\ 

at every point on the grid. Here r is the radius, v is 
the fluid velocity vector and | is the absolute value of 
velocity^. As the shock strongly compresses the g as, a 



Fig. 3. — Post outburst: the radial profile of the physical vari- 
ables when a burst appears at Ati = 0.235 Myr. Solid, dotted 
and dashed lines correspond to At = 0.235, 0.239 and 0.243 Myr. 
The fiuid radial velocity changes direction from inward to outward 
between Ati = 0.235 Myr to At2 = 0.239 Myr. The cold shell is 
still visible in the solid line at the center. The arrows (A, B, C) in- 
dicate the positions of the shock at the three defined moments (see 
Fig. |4jl; the shock initially accelerates (A— >B) due to the SMBH 
driving from below and then decelerates (B— >C) as adiabatic losses 
overcome the SMBH input. 

local deep minimum of S marks the position of the shock. 
The upstream sonic Mach number of the shock is given 

by 



Cs 

where Us is the shock velocity relative to the galaxy and 
Cs is the upstream sound speed. The shock velocity Us is 
numerically estimated by determining the change in the 
shock position within a time step. The upstream sound 
speed is calculated from the upstream temperature T and 
density p as 




where P is the pressure, ks is the Boltzmann constant, 
fi = 0.62 is the mean molecular weight, rup is the proton 
mass, and 7 = 5/3 is the adiabatic index of the gas. 

As the magnetic field B, necessary for the synchrotron 
emission, is not included in our hydrodynamical code, 
we must assume a fiducial value based on physical ar- 
guments. In particular, for elliptical galaxies, there is 
evidence that the ratio between the n on-thermal and 
therm al pressures is about 0.1 — 0.2 (e.g.. lChurazov etldl 
I2OO8I : Churazov et al. in preparation, private commu- 
nication). Here we just assume the ratio between the 



^ Only the radial component is non-zero due to adopted spherical 



symmetry. 
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thermal and magnetic pressure to be 

P 



10, 



(4) 



so that we can estimate the postshock magnetic field from 
the thermal pressure. 

Now we summarize the basic results of the test par- 
ticle theory. The 1-D momentum spectrum of the non- 
thermal electrons is taken to be a power law : 



f{p) = 4V X Co P = 47rCo p 



2-q 



(5) 



where p is the momentum of the electrons and Cq is the 
normalization of the spectrum, which will be determined 
later. The exponent q is given by^ 



3t 



1 ' 



(6) 



where the compression ratio r is related to the Mach 
number M as 



7- 



1 



1 



7-1-1 7-1-1 



(7) 



From the momentum distribution ([5]), and from the 
energy-momentum relation 



p = m^c X 



1 



E 



-1, 



(8) 



where E and me are the kinetic energy and the rest mass 
of electrons respectively, and c is the speed of light, we 
can express the energy spectrum of the non-thermal elec- 
trons as 



N„{E)^4ttCo p 



dE 



(9) 



The energy of non-thermal electrons ranges between 
Emin (with the corresponding minimum momentum 
Prnin) and E,nax (with maximum momentum p,nax)- 
Following what is done in t he computation s of SNRs 
and galaxy clust ers (e.g., [Berezhko & EUisonl 119991 : 
lEnBlin et al.l[200l . the minimum momentum Pmin is re- 
lated to the downstream temperature T2 by the injection 
parameter Xmj'- 



Pn 



^171 J 



\/2xm^kBT2 (cgs). 



(10) 



For electrons, the maximum energy at the shock is deter- 
mined by the balance betwee n radiative loss and accelera- 
tion g ain, which gives (e.g. , IWebb et al.l[l98l IRevnoldd 
I2OO8 



E„ 



= 71.65 (w,/10**) X (Ba/lO"")""-^ erg (cgs) 



= 44.72 7/,,io3k„Vs X (^2,^0)"*^" TeV 



(11) 



^ Here we neglect the Alfven speed, with which the scattering 
centers presumably move, in comparison to the flow speed. This 
correction is only important for low Mach number shocks. Because 
synchrotron emission mainly comes from strong shocks, ncglcct- 
ing Alfven sp eed will not change the spectrum significantly (e.g., 
ICaprioh et~al. 2009). 

* In a typical elliptical galaxy, the energy density of IR back- 
ground and star light is usually small compared with the down- 
stream magnetic energy density. So we ignore the electrons' inverse 
Compton cooling here. 



where B2 is the downstream magnetic field. Actually, 
in the test particle theory, the synchrotron spectrum is 
not very sensitive to the maximum energy Emax- But 
the cutoff of the power-law spectrum does depend on 



En 



In order to determine the normalization of the 



electron spectrum, we follow the a pproach adopted in 
simulations of galaxy clusters (e.g., lEnBlin et al.l l2007t 
iPfrommer et al.ll2008f ) and smoothly connect the power- 
law momentum spectrum of electrons to the downstream 
thermal spectrum at the momentum Pmin- Then the 
coefficient Co in equation ([5|) is given by 



Cn = 



(2^mefcBT2)3/2 



-P™i„/(2mefcf3 7^2) 



(12) 



where rig, 2 is the downstream electron number density. 

The test particle result is valid only when the energy 
density of the non-thermal particles is small enough com- 
pared to the thermal energy density so that the CRs do 
not significantly change the structure of the shock. How- 
ever, for a strong shock with a very high Mach number 
(e.g., when the shock encounters the cold shell), calcula- 
tion with the above method will easily produce an elec- 
tron energy density comparable to the downstream ther- 
mal energy. So in our calculations we have to choose ap- 
propriate normalization of the electron momentum spec- 
trum for differ ent Mach numbers , an d we adopt the 
metho d used in lEnBlin et al.l (|2007| ) and IPfrommer et al.l 
(|2008h . First, we obtain the electron energy spectrum 
from equations ^ and (fT2|) . so that the energy density 
of the non-thermal electrons is 



ENo{E)dE . 



(13) 



Then, the ratio between the non-thermal electron energy 
density and the downstream thermal energy density (in 
the test particle theory approach) is 



Clin 



K, 



e,0 



1.5n2kBT2 



(14) 



where n2 is the downstream total number density. The 
ratio £,iin can be very large for a shock with large M. 
Now we set a maximum ratio between the non-thermal 
electron energy density and the downstream thermal en- 
ergy densi ty to £.max =0-05 ( a value s uggested by sim- 
ulatio n of iKeshet et al.l 120031 : also see IPfrommer et al.l 
l2008f ) and define a parameter x = iun/imax- The mod- 
ified normalization of the electron momentum spectrum 
we use in our calculation becomes 



C = (l-e-x)x-^Co 



(15) 



so that C ^ Co for x ^ 0, while C Cq/x f^or very 
large values of x- In other words, if S^un is very small 
compared with S^max (x ^ 1); then C ~ Co and the energy 
density ratio is approximately S^un- However, if is 
very large compared with £,rnax (x ^ 1)j the final energy 
density ratio between the non-thermal electrons and the 
downstream thermal energy density cannot exceed £,max- 
Roughly speaking, the modification is important when 
M is larger than 10. Therefore, equation © becomes 



N{E) 



2-q dp 



(16) 
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and the energy density and number density of the non- 
thermal electrons are given by 



K,. 



EN{E)dE ,iVe 



N{E)dE.{n) 



The finally ratio between the non-thermal electron en- 
ergy density and downstream thermal energy density ^ 
is the same as equation (fT4|) with Ke,Q replaced by Ke- 
In particular, the injection coefficient 77, a parameter de- 
termining the acceleration efficiency model ed in the dif- 
fusive shock acceleration mechanism (e.g., lEnBlin et al.l 
I2007f ). is defined here to be 



ne,2 



(18) 



Another thing which must be emphasized is that here 
we do not give the energy density of the non-thermal 
protons, which is unnecessary for our purpose. By us- 
ing test-particle theory we are assuming non-modified 
shocks, and the CR energy density is less than 10% of 
thermal energy density by assumption^. This is only a 
first order approximation and a first step to determine 
whether a more advanced model is merited. 

2.3. The Spectrum after the Shocks 

The test particle framework described above gives the 
energy spectrum of the non-thermal electrons at the 
shock, and we need to compute the evolution of the spec- 
trum further downstream. For simplicity, we assume no 
diffusion in our calculation. This means that once the ac- 
celerated electrons are generated at the shock, they are 
frozen to the fluid element. So, while the shock propa- 
gates outwards, the non-thermal electrons are left behind 
and move with the fluid. These non-thermal electrons 
suffer radiative and adiabatic losses, and the electron en- 
ergy spectrum will change with time. The changes of 
the el ectron energy and number density are well- known 
(e.g.. IRevnoidi [l99l [Cassam-Chenai et al.l [20071 ). We 
will describe our numeric al implementatio n of the spec- 
tral evolution formulae of ' Revnoldsl (|1998l ). We assume 
that an electron with energy Eq is generated at time t^ 
at radius tq, where the mass density is po- At a later 
time the electron moves to the new position r, where 
the mass density is pt = p{r, t) , so that the compression 
ratio of the fluid element is 



a{r, t) = — . 
Po 



(19) 



After the losses, the energy of the electron at time t and 
radius r is 



E{r,t) = a 



= ^1/3. 



En 



l + BEo 



(20) 



where the quantity Q is the result of the following inte- 
gral 

Q{t)^Af Bl^{t)a^/\t)dt. (21) 

Jta 

^ Some observations show that in SNRs 10% of the mechan- 
ical energy of the explosion can be converted into CRs (e.g., 
rSiaronian et al..,2004i') . 



The coefficient in front of the integral is ^ = 
4e'^/(9mgC^) — 1.57 x 10^^ (cgs) (e is the electron's 
charge ) and 



(22) 



where i?cbr = 3.27 is the magnetic field strength with 
energy density equal to that of the CMB at the redshift 

2 = 0. 

Each of the non-thermal electrons generated by the 
shock will change its initial energy according to the above 
equations, and the energy spectrum of these electrons 
when they are generated is given by the test particle 
result (equation [TB)) . As the number of electrons is con- 
served, at time t and radius r the energy spectrum of 
these electrons becomes 



N{E)r^t = NiEo)a 



-2/3 



Eo{E) 



E 



(23) 



Note that Eq as a function of E is determined by in- 
verting equation (j20|l . In particular, initially the max- 
imum energy and minimum energy of the electrons are 
Emax.o and Emin,o- Subsequently, the maximum energy 
and minimum energy will also change to the new maxi- 
mum energy Emax and minimum energy Emm according 
to equation ([^ . 

Once we obtain the energy spectrum of the electrons, 
wc can calculate the resulting synchrotron spectrum, i.e., 
the emissivity of the population of electrons in the energy 
interval Emin < E < Emax 



Viiy) N{E)dE erg s~icm"3Hz^^(24) 



Then the luminosity L^, can be calculated by integration 
of Jjy over the whole volume. In the above equation, the 
synchrotron emission from a single electron is given by 



K5/z{x)dx erg s"^Hz"^(25) 



c J v/vc 



where 1^5/3 (a;) is a modified Bessel function of the second 
kind, and 



47rm2c^ 



(26) 



is the critical frequency, where E is the energy of the 
electron and B±^ is the magnetic field component per- 
pendicular to the lin e of sig ht (e.g.. lBlumenthal fc Gouldl 
1970; Lazendic et al.ll200l . In our calculation B± is ap- 
proximated by B given in equation ([4]) . 

2.4. Numerical Method 

As the shock emerges from the center of our model 
elliptical galaxy and moves outwards, it will accelerate 
particles at each position it passes through. Once the 
non-thermal particles are generated by the shock, they 
will be frozen in the fluid element as we assume no diffu- 
sion, and the spectrum will evolve according to equation 
(f23|) . In practice, in our modeling, at every given time, 
there will be non-thermal electrons within a certain ra- 
dial range with the new ones produced by the shock and 
the old ones moving out with the fluid. Accordingly, we 
calculate the synchrotron emissivity from equation l|24p 
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at each radius and then obtain the total luminosity by 
integrating over the volume. 

In practice, once the electrons are generated by a shock 
at certain radius tq at time Iq, we follow the fluid element 
and compute its position at time t. From the hydrody- 
namical code, we know the velocity Vr{t) and acceleration 
ttrit) of each fluid element, so that the new position rt 
can be estimated as 



6xl0» 



(27) 



where the position increment Ar in each time step At is 
calculated as 

A r = Vr{t) X At+ ^ar{t) x At"^ , (28) 

For higher numerical accuracy, the quantities Vr{t) and 
ar(t) of the fluid element at radius r could be approxi- 
mated by the average velocity and acceleration of the two 
fluid elements at radius r and r + dr, which are given in 
the output from the code. 

So, in our calculation, when a shock is found at time t, 
its position is determined by the criterion given in equa- 
tion ll]). Then we follow the fluid element, that at time 
t was at the position of the shock according to equa- 
tions (P7|) and After one time step, we identify 
the new position of the shock and follow a new fluid el- 
ement. While following each fluid element that passed 
the shock, we calculate the value of Q(t) given by equa- 
tion (j2ip . In this way, we know the energy spectrum 
of the electrons from equation ([23]) at different positions 
at time t, and we flnally calculate the total synchrotron 
spectrum. Note that our numerical method is different 
from the standard method which solves the continuity 
equation for el ectrons in the 2-D phase space at each 
time step fe.g.. iLongaiii [l99^ . Instead, in our method, 
the energy spectrum at each time step is calculated ana- 
lytically according to equation (|23|) , and we only need to 
follow one-dimensional the fluid element, which is com- 
putationally more efflcient. 

3. RESULTS FROM ONE SHOCK 

Now we give an example of the determination of the 
properties of the synchrotron spectrum for a shock in 
model B3o2. In particular, we focus on the last major 
burst in the example described in ^2.1\ 

3.1. Synchrotron Spectrum 

We define the time when the SMBH bolometric ac- 
cretion luminosity Lbh reaches 1 percent of Eddington 
luminosity LEdd as the zero time point, corresponding to 
6.48 Gyr in the code. Before the shock is formed and 
can be recognized by the code, a cold shell is falling to- 
wards the center, which can be seen clearly from Fig. [21 
From Fig. [31 we can see that at a later time At = 0.235 
Myr, the radial velocity changes from negative to posi- 
tive, indicating the emergence of the nuclear wind and 
the shock. The result of equation ([T]) is shown in Fig. 
[H At the position of the shock, the parameter S has a 
minimum value, which confirms our criterion for finding 
the shock. 

At each time we find the position of the shock and 
follow the shock until it disappears. Meanwhile, we 
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Fig. 4. — The shock-locating method. In the two top panels, the 
dimensional quantity div(v) (left) and the dimensionless function 
5 (right, equation [TJ are shown for the same profiles shown in Fig. 
\3\ The solid, dotted and dashed lines are in time order, separated 
by 0.004 Myr. The solid line is at time Ati = 0.235 Myr after the 
time zero point. The three lines also show the movement of the 
shock, which can be used to estimate the shock velocity. There are 
three arrows and letters in the figure to show the position of shock 
found at different times. 



calculate the energy spectrum of the non-thermal elec- 
trons accelerated by the shock with equation ([TB|) and 
the spectral evolution according to equation Then 
the synchrotron emissivity at each radius is calculated 
from equation (j24[) and the total synchrotron spectrum 
for our model elliptical galaxy is the sum of the emission 
from different positions. So, the final spectrum is actu- 
ally not from a single population of power law electrons, 
but from electrons at different times in their evolution. 
In Fig. show the radial profile of the magnetic field, 
the non-thermal electron energy density and the flux for 
different frequencies at the time At — 0.335 Myr. The 
magnetic field is estimated according to equation ^ . At 
this time, the shock is located at about 327 pc and the 
bump of the magnetic field is due to the compression at 
the shock. The non-thermal electron energy density is 
estimated according to equation ([17]). As we assume no 
diffusion, non-thermal electrons only exist in the radial 
region 258 — 327 pc at this time. The non-thermal elec- 
trons around 260 pc come from the initial strong shock 
(arrow A in Fig. [J). The flux at radius R is estimated 
by vJi, X R, and the fluxes of radio emission (4.84 GHz), 
optical emission (7.25 x lO^"' GHz) and X-ray emission 
(2.42 X 10^ GHz) are labeled by sohd, dashed and dotted 
lines, respectively. Though the electron energy density 
around 260 pc is large, these are "old" electrons that do 
not contribute to the synchrotron emission in the three 
frequencies by this time. The X-ray emission can only be 
found just behind the shock, while the optical emission 
can be found within two cells of the shock. However, the 
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Fig. 5. — Radial profile of tlie magnetic field, electron energy density and synchrotron flux for three frequencies at time Ai = 0.335 
Myr. At this time, the shock is located at about 327 pc. The upper panel shows the radial profile of magnetic field estimated according 
to equation (|4jl at this time. The bump around 300 pc is d ue t o the compression by the shock. The middle panel shows the radial profile 
of electron energy density calculated according to equation l|17|l . As we assume no diffusion, non-thermal electrons only exist in the radial 
region 258 — 327 pc at this time. The lower panel shows the radial profile of the fiux for different frequencies: radio band 4.84 GHz (solid 
line), optical band 7.25 X 10^* Hz (dash hne), X-ray band 2.42 X 10^'' Hz (the dot). The flux at radius R is estimated by uJv X R. Note 
that the higher the frequency is, the closer the fiux is to the shock. At this time, the X-ray emission can only be found just after the shock, 
while the optical emission can be found in two cells after the shock. For the radio emission, it exists in a much more extended region. 



radio emission can extend over almost the whole region 
passed by the shock. 

In our model, the only free parameters needed to cal- 
culate the synchrotron emission are Xmj , defined in equa- 
tion (fTU|) . which is related to the minimum momentum 
of the accelerated electrons, and the saturation param- 
eter ^max = 0.05. In order to make the injection co- 
efficient rj (equation [T8| to be about lO""*, which is 
a reasonable value used in the models of SNRs (e.g., 
iBerezhko &: Ellisoiilll999f ). we take the valu e^ Xjni = 3.6, 
the same as the value taken in the model of lEnfilin et al.l 
([2007). Note that as we do not include the jet feedback 
in the hydro simulation (see Paper I), the synchrotron 
emission calculated here is not the total emission that 
would be observed from a real elliptical galaxy. It can 

^inj is smaller, more electrons will be accelerated and the 
emission will be stronger while the contrary is also true. However, 
as electrons with momentum around Pmin do not contain most of 
the energy for non-thermal electrons in strong shocks, the total 
spectrum is not very sensitive to Xi„j. 



only be compared to elliptical galaxies without jets or 
to the core emission in elliptical galaxies with extended 
jets. We take snapshots at four different times during 
the evolution of the model. The resulting spectra are 
shown in Fig. [SI The change of spectral shape with time 
in this figure is due to the radiative and adiabatic losses. 
In particular, the high frequency part disappears at late 
times (due to a smaller magnetic field at larger radii and 
losses). We define the fiducial radius rh enclosing the vol- 
ume within which half the energy at the frequency 5 GHz 
is emitted. The spectral index a at a certain frequency is 
defined to be the index of the power law approximation: 

J^(xv" . (29) 

In practice, we calculate the spectral index by the fol- 
lowing difference 

log J^, - log J^, 
a = -j , (30) 

log 1^1 - log 1^2 

where the two frequencies are chosen to be vi — 0.17 
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Fig. 6. — Synchrotron emission spectra without absorption at four different times. The soUd, dash, dot and dash-dot lines are the total 
synchrotron spectrum from the elliptical galaxy at the time Aii = 0.235 Myr, At2 = 0.275 Myr, Ata = 1.035 Myr and At4 = 4.635 Myr 
respectively after the defined time zero. The half radius w ith each line is the radius within which half energy at the frequency 5 GHz 
is emitted, a is the spectral index defined in equation II30I I. The arrow labeled with 'u = 5 GHz' is the range of observed core radio 
luminosity at 5 GHz in FR H sources llZirbel Sz Baumlll995l : [Chiaberge et al.ll2000l) . The arrow labeled with = 4.3 X 10^* Hz' is the range 
of observed core optical lu minosities at 4.3 X 10^^ Hz (wavelength A = 7000A) in FR II sources by HST, which are believed to originate 
from non-thermal emission llChiaberee et ani2000l) . At the bottom of this figure, the names for different frequency bands are labeled. 



GHz and 1^2 = ^ GHz, so that we can compare our calcu- 
lation with observations in the two bands. From Figure 
[6l we can see that the spectrum becomes steeper with 
time and the absolute value of the spectral index be- 
comes larger. The ratio ^ of the non-thermal electron en- 
ergy density generated by the shock to the downstream 
thermal energy density will change for different Mach 
numbers, as shown in Fig. [71 The non-thermal electron 
energy density is calculated according to equation p?]) . 
Note that the maximum ratio is ^ ~ 0.05, which is nearly 
reached only for very high Mach numbers. For small 
Mach numbers, the ratio is only about 0.001. Roughly 
speaking, when M is larger than 10, the saturation con- 
dition (equation llSp becomes very important. Though 
we do not consider protons here, they will also be accel- 
erated by the shock. As total cosmic ray energy density is 
estimated to be around 10 percent o f the shock's kinetic 
energy (e.g.. lAharonian et al.ll200^ . the energy density 
ratio between the protons and electrons will also change 
with Mach number. For large Mach number, the en- 



ergy density ratio between the non-thermal protons and 
non-thermal electrons is smaller (say 10 or smaller) while 
for small Mach number, the ratio is larger (for example 
100). The adopted test particle theory assumes that the 
post shock CR pressure is small compared to the ther- 
mal gas pressure and we find that this assumption to be 
valid during most of the time of evolution except when 
the shock encounters the cold shell, where the modified 
spectrum takes effect. To the first order of approxima- 
tion, we neglect the pressure of CRs, which is a first step 
to determine whether a more advanced model is merited. 

3.2. Synchrotron Radio Emission 

While the shock moves outwards, it dissipates energy 
into the ISM and finally disappears. In this example, we 
can still see the shock at the time At — 46 Myr at the 
radius 32.7 kpc, where a second shock due to the next 
sub-burst emerges. But as we have shown in Fig. [SI the 
radio emission is already almost unobservable at the time 
At = 4.635 Myr. In Fig. [51 we show the evolution of the 
synchrotron emission size r^ (at 5 GHz, defined as the 
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Fig. 7. — Acceleration efficiency at different Mach numbers, A4. 
The figure shows the ratio between the non-thermal electron en- 
ergy density accelerated by the shock to the downstream thermal 
energy density just after the shock ^, as computed from the hydro- 
dynamical code. For high J\4, the ratio is almost 0.05, which occurs 
when the shock encounters the cold shell. For small Ai, the ratio 
can be around 0.001. To the degree that ^ <S; 1, the test particle 
theory adopted in this paper is accurate. Roughly speaking, when 
A4 is larger than 10, modification to the normalization (equation 
list is important. The red dashed line is the theoretical result of the 
model with an assumed constant downstream temperature 5 X lO'^ 
K, an assumed constant downstream number density 10 cm~^, and 
an assumed maximum electron energy 100 TeV. Even for the same 
A4 the ratio ^ can be different for different temperature and num- 
ber density. 
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Fig. 8. — Evolution of the synchrotron radio emission size at 5 
GHz. The size Tf^ is defined to be the radius within which half the 
energy at the frequency 5 GHz is emitted. At around 3 Myr, the 
size decreases because at this time there is no contribution to the 5 
GHz synchrotron radio emission from the new shock position due 
to small magnetic field. Then the emission at this band is only 
from the 'old' electrons, which suffer large radiation and adiabatic 
losses and move outwards. 

spatial radius enclosing half of the emitted synchrotron 
luminosity at the given time). Generally, rh increases 
with time as the non-thermal electrons, which are frozen 
in the fluid element, move outwards. However, when the 
shock propagates into the regions with low magnetic field 
strength, there will be almost no synchrotron radio emis- 
sion at 5 GHz and the apparent size will decrease at that 
time. The size rh may change from a very small initial 
size (5.1 pc in this example) to around 4.5 kpc. In Fig. 
[9l we show the evolution of synchrotron luminosity at 3 
different frequencies. Note that the spectrum shown in 




Fig. 9. — Evolution of the synchrotron luminosity in 4 different 
bands: v = 0.17 GHz (solid dots, radio), u = 4.84 GHz (crosses, 
radio), u = 7.25 X 10^ GHz (triangle, optical, and the first two 
points in this line are connected by a dash line because they will 
be absorbed by the dust), u = 2.42 X 10^ GHz (open circles, X- 
Ray). The larger the frequency is, the sooner radiation at that 
frequency disappears. In the four horizontal dashed lines we show 
what would be seen if uLi, for each cited frequency band had the 
power 5 X 10^* erg s~^. 

Fig. [6] is the emission spectrum without allowance for 
absorption. However, for radio bands, the optical depth 
due to synchrotron self-absorption at 0.1 GHz initially 
is only about 10~^ and smaller at later times. The syn- 
chrotron self-absorption is thus only important for fre- 
quencies smaller than 10^ Hz in this example. 

3.3. Synchrotron Emission in Other Bands 

For optical bands, the main absorption is due to dust. 
The optical depth due to dust absorption is calculated 
in the code (Paper I). The optical emission is absorbed 
by dust for At < 0.33 Myr (within 300 pc in this ex- 
ample), but the absorption is not important after that 
time. A fiducial estimate of the expected mean flux 
Eo of the optical synchrotron emission is obtained as 
Eo = Lo/(27rr^g), where Lo is the total synchrotron 

luminosity in the optical band (7.25 x lO'^'' Hz), and 
where we have neglected the difference between the ra- 
dius within which half of luminosity is emitted r^^o and 
the projected surface effective radius. We compare Eq 
with the optical flux of our galaxy model, to see whether 
the synchrotron optical emission can be observed or not. 
In order to estimate the stellar flux of the galaxy model, 
we adopt as a flducial value for the surface brightness 
at the effective radius Eg = 22 mag/arcsec^, and for the 
effective radius Re = 6.9 kpc. At each time, we then com- 
pare the value of Eq with the stellar flux profile obtained 
from the de Vaucouleurs profile, i.e., with E^(r;i_o/i?e). 
The result is shown in Fig. [TUl This figure shows that 
the synchrotron optical emission is stronger than the star 
light only at early times. Later on, it is too faint com- 
pared with the star light. However, the optical emission 
will be absorbed by the dust in this example at the early 
time. So we expect to observe synchrotron optical emis- 
sion at the inner part of those elliptical galaxies with little 
dust during a short time in the quasar phase. This result 
is just from one shock at the first sub-burst of the last ma- 
jor burst in Fig. [TJ However, we stress that within each 
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Fig. 10. — Properties of synchrotron (optical) emission in model B3o2. This figure shows the time evolution of the flux S (upper panel), 
luminosity Li, (middle panel) and size „ (lower panel) at the frequency 7.25 X 10^"* Hz. The dashed line in the top panel is 12i,{rh,o/ Re), 
i.e., the stellar flux of a de Vaucouleurs profile with = 22 mag/arcsec^ and Rf, = 6.9 kpc. In the upper and middle panels, the first two 
points are connected by dotted lines because the optical emission will be absorbed by dust at that time. Note that synchrotron emission 
lasts for 2.8 Myr and then drops to zero very quickly. 



burst several shocks appear from the sub-bursts in the 
galaxy almost simultaneously, which will give stronger 
emission and a higher probability of observing the syn- 
chrotron optical emission. The study of the cumulative 
effect of several co-existing shocks is beyond the reach of 
the present investigation, and it is postponed to a future 
work. In particular, it will be important also to con- 
sider the more realistic case of multiple shocks created by 
the so-called "time-dependent nuclear wind treatment" 
which includes the time for the wind to reach different 
radii (sec Paper I , see also Sj4]). Right now, we have just 
checked that, for a single shock, the time-dependent and 
time-independent wind (neglecting the time needed for 
the wind to reach different radii) models do produce very 
similar results. 

Fig. [6] shows that there is also synchrotron X-ray emis- 
sion for the duration of ~ 10"* yr to ~ 10^ yr at the 
very early time of the burst. At the time At — 0.235 
Myr, at the frequency 2.4 x 10^^ Hz (equivalent 1 keV), 
the luminosity is 5.76 x 10^^ erg while at the fre- 
quency 1.2 X 10^^ Hz (equivalent 5 keV), the luminosity 
is 1.2 X 10^^ erg s~^. The synchrotron X-ray emission 



lasts for only 0.32 Myr and drops to zero rapidly af- 
ter that time. Typically the luminosity of the thermal 
X-ray emission in the band 0.5 — 2.0 keV from isolated 
elliptical galaxies va ries from 10"^^ to 10"*^ erg s~^ (e.g., 
iMemola et"all l200tf). So the synchrotron X -ray emission 
are likely to be observed only in these "on" state galax- 
ies, which are only a small fraction of all galaxies. Be- 
cause the spectrum of bremsstrahlung emission (energy 
density per Hz) is almost flat, we should expect the lumi- 
nosity of bremsstrahlung optical emission to be around 
10"^^ to 10"^^ erg s~^ based on the bremsstrahlung X-ray 
emission. This is actually smaller than the synchrotron 
optical emission here. 

As synchrotron X-ray emission is also observed from 
SNRs, we can scale down the flux of the X-rays calcu- 
lated here to a typical flux of X-rays from SNRs. If syn- 
chrotron emission from SNRs and elliptical galaxies is 
from the same mechanism and we keep the ratio of syn- 
chrotron optical emission and synchrotron X-ray emis- 
sion calculated here, we can estimate the flux of syn- 
chrotron optical emission from SNRs and see if it can 
be observed. Because the synchrotron spectrum basi- 
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Fig. 11. — Properties of the synchrotron (soft X-ray) emission in model B3o2- This figure shows the time evolution of the flux S (upper 
panel), luminosity (middle panel) and size r^ ,^ (lower panel) at the frequency 2.41 X 10^^ Hz at different times. The half radius in the 
X-ray band is almost the same as the radius of the shock, which means the synchrotron X-ray emission only exists within a narrow region 
around the shock. 



cally becomes steeper with time, we can use the ra- 
tio at the initial time to estimate the lower bound of 
the synchrotron optical emission from SNRs. The syn- 
chrotron X-ray flux obser ved from SNR G 330. 2-1-1.0 is 
about 10"" erg s^^ cm'^ (jPark et al.ll2009l ). In our cal- 
culation, at the time At = 0.235 Myr (see the solid line 
in Fig. O, the X-ray flux (at 1 keV) is about 1.85 x 10^ 
erg cm~^ while the optical flux is about 2.83 x 10^ erg 
s~^ cm~^. Then we give a lower bound of synchrotron 
optical emission from SNRs to be 1.53 x 10~^^ erg s~^ 
cm~^, corresponding to 44.87 mag/arc^ without absorp- 
tion. Actually, in SNR G332.5-5.6, line emission with a 
flux as faint as 10~ ^^ erg s~^ cm~^ has been observed 
(jStupar et al.l 120071 ). It seems that the synchrotron op- 
tical emission from SNRs should be observable by the 
appropriate technique. 

4. DISCUSSION 

4.1. Observational Consequences of Our Model 

Our model has observational predictions, some of 
which are consistent with current radio observations and 
some of which may be feasible Chandra source detec- 
tions. First, as shown in Fig. [6l the synchrotron emis- 



sion at 5 GHz and 4.3 x lO^"' Hz falls in the range 
of observed values in FR II sources, which means that 
our mechanism should not be neglected. Besides, re- 
cent observations at 1.4 GHz find some early- type core 
galaxie s hosting radio-loud A GN of extremely low radio 
power (|Baldi fc CaDetti|[2009l ). which can be interpreted 
as radio emission via the mechanism described here as it 
would be too weak to be produced by typical jets. Sec- 
ond, as shown in Fig. [8l the shocks produced in the sim- 
ulations easily reach the kpc scale before they die away 
gradually. The synchrotron radio emission can also go to 
this large scale (though it is very weak). Actually, kpc- 
scale outflows or w inds have been obse rved in hig h red- 
shift quasars (e .g., !de Kool et al.l 120011: iNesv adba et aT] 
120081: INesvadba 2009; Cattaneo et al.ll2009f ). which pro- 
vides strong support for our scenario based on the com- 
puted winds in th e AGN simulations. I n particular, a 
recent observation (Ale xander et al.ll2009( ) finds a galaxy 
wide outflow at redshift ~ 2, which is believed to be 
a wind radiatively driven by the AGN and/or super- 
novae winds rather than by jets. N ote that some ob- 
served X-shaped radio galaxies (e.g., Leahy fc WilHamsl 

" 20081 ) provide sig- 
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nificant support for the idea that the observed radio 
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emission arises from shocks caused by feedback from the 
SMBHs, as the broad Hne winds drive shocks into the 
ambient gas with a similar s hape (e.g., iDrew fc Progal 
[2OOOI iDorodnitsvn et al.][200l l. 

Third, the duration of synchrotron emission at a cer- 
tain frequency will decrease with increasing frequency. 
For a single shock, synchrotron X-ray emission can only 
exist behind the shock, where the non-thermal electrons 
are just created. The synchrotron optical emission can 
extend to a larger region while the synchrotron radio 
emission can almost extend to the whole region where 
the shock has passed, which can be seen from Fig. [5l 
So, in principle, at good enough resolution we should see 
the synchrotron emission in different bands at different 
spatial positions. 

Though synchrotron emission from our mechanism and 
jets is similar to a certain extent, there are some differ- 
ences that can help distinguish the two mechanisms. In 
the following, we list four possible observational tests. 
First, we expect the magnetic field in the shocked ma- 
terial to be more randomized than the magnetic field in 
the jets. So, the polarization of the emission from our 
mechanism will be weaker than the polarization of the 
emission from the jets, and polarization observations can 
help to disentangle the two mechanisms. Second, as the 
shock (and the synchrotron emission) is produced after 
each AGN outburst, emission from our mechanism will 
only exist in a certain fraction of elliptical galaxies, pro- 
portional to the duty cycle. In practice, in a survey of 
elliptical galaxies, where the emission from jets gives the 
background level of core radio emission, the number of 
galaxies with extra radio emission should be comparable 
to the duty cycle (see Paper I for a discussion of duty 
cycle of the models). Third, as the velocity of the jet is 
typically much larger than the velocity of the shock, then 
if we can directly measure the velocity of the outflow, it 
will be easy to tell which process dominates. Finally, 
the variability time scale of emission from jets is shorter 
than that of our mechanism. The variability time scale 
can also help us tell the origin of the core radio emission. 

4.2. Possible Future Improvements of the Model 

There are still some significant uncertainties in our 
model and further improvements are needed. First, the 
hydrodynamical code does not include the magnetic field, 
a fundamental ingredient for the accurate computation 
of synchrotron emission. In this paper, we only assumed 
a ratio between the magnetic pressure and thermal pres- 
sure to estimate the preshock magnetic field in elliptical 
galaxies. A better method is needed to give the distri- 
bution of magnetic field and its evolution with time. 

In fact, the magnetic field will be amplified near the 
shock due to back reaction of cosmic rays, and we roughly 
estimate that the amplified magnetic pressure should 
not exceed about 10% of the postshock thermal pres- 
sure. A mo re careful calculatio n of this effect (e.g . . iBel] 
2004 iRiau elme & Spitkovsky 2001 lOhira et all 120091 : 
Zirakashvih fc Ptuskin II2008I) would be preferable. Sec- 
ond, we use the test-particle theory to calculate the en- 
ergy spectrum of the non-thermal electrons and neglect 
the modification of the shock structure due to the cosmic 
rays. However, as is already known in SNRs, non-linear 
theory is needed to calculate the spectrum of the non- 
thermal electrons for strong shocks, especially at high 



energ ie s fe.g.. iBerezhko fc Ellisonlll999l: [Lazendic et al.l 
2OOI iLee et al.l l2008t iTreumann fc Jaroschekl 120081 : 
Revnoldsll2008ir So in order to calculate the synchrotron 
emission more accurately, we will need to include the 
diffusive shock acceleration mechanism into the hydrody- 
namical code. This greatly adds to the complexity of the 
modeling. Third, as anticipated in ii2.11 all the results 
presented here are obtained using a 1-D code. Of course, 
this is quite a strong approximation, which precludes the 
quantitative analysis of several phenomena that almost 
certainly take place in real galaxies. One of the most im- 
portant in the present context is the possible cold shell 
fragmentation due to Rayleigh- Taylor instability follow- 
ing central outburst. While we could not exclude this 
possibility, we note that several properties of our 1-D 
models (such as the duty-cycle, mass accreted on the 
central SMBH, star formation. X-ray luminosity) agree 
well with the observations, so that we are confident that 
the basic physics of feedback is actually captured by our 
simulations. Note also that on the kpc scale, the cen- 
tral outbursts, even if initially anisotropic, will necessar- 
ily become rounder and rounder as they propagate, and 
therefore better and better approximated by our code. 

The calculations here are based on the so-called "time- 
independent nuclear wind model" (with mechanical as 
well radiation feedback). In this simplified approach, 
which saves much simulation time, the propagation ve- 
locity of the wind along the numerical grid (Paper I, eq. 
29) is set to be infinity. Therefore, the nuclear wind mass, 
momentum and kinetic energy are instantaneously dis- 
charged over the whole computational grid. The more re- 
alistic model is the time-dependent nuclear wind model, 
which will be more accurate, especially at the inner part 
of the galaxy at early times (§2.3 of Paper I). More ac- 
curate calculations of the synchrotron optical and X-ray 
emission also require the time-dependent wind model. 
Finally, effects of multiple shocks on the emission from 
the inner galaxy and on the reacceleration of preexisting 
nonthermal particles have yet to be included. 

These caveats aside, the presented results show that 
central outbursts of mechanical and radiative energy con- 
sequent to AGN flaring should produce shock accelerated 
electrons capable of producing the non-thermal emission 
actually observed from the cores of some elliptical galax- 
ies subsequent and consequent to AGN outbursts. 

5. CONCLUSIONS 

In this paper, we have explored the observational con- 
sequences of a new mechanism for producing the syn- 
chrotron emission seen in elliptical galaxies. This idea 
is based on particle acceleration subsequent to bursts of 
energy from the central black holes. The properties of 
the ISM are computed from hydrodynamical simulations 
including radiative and mechanical feedback consequent 
to an AGN outburst. Due to this feedback, a wind is 
formed during the quasar phases, and a shock is driven 
into the elliptical galaxy. Given standard physics, the 
shock, during its propagation into the galaxy, will accel- 
erate electrons and protons via the diffusive shock ac- 
celeration mechanism. Synchrotron emission arises from 
these non-thermal electrons. In this paper, we have fo- 
cused on the computation of the energy spectrum of the 
integrated synchrotron radiation. We include radiation 
loss and adiabatic loss during the evolution of these non- 
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thermal electrons and we show the change of the syn- 
chrotron spectrum with time. 

In general, we have found that the synchrotron radio 
emission from one shock can last about 4.7 Myr at 5 GHz, 
while a shock can exist for 46 Myr. During the evolution 
of the elliptical galaxy, there are several major bursts and 
several sub-bursts for each major burst. Roughly speak- 
ing, the shock can exist during the whole burst phase of 
the galaxy while the radio emission (around 5 GHz) ex- 
ists for 10 percent of the burst phase. Besides, our mech- 
anism also gives marginally detectable synchrotron opti- 
cal emission during the 3 Myr after a burst. The duration 
of synchrotron optical emission is about 50% of the radio 
emission at ~ 0.17 GHz while the duration of the syn- 
chrotron X-ray emission is approximately 5% of the radio 
synchrotron emission at ~ 0.17 GHz. The half radius of 
the synchrotron X-ray emission is almost the same as the 
radius of the shock, which means that the synchrotron 
X-ray emission only exists as a 'ring' around the shock 
while the half radii of the synchrotron radio emission 
and synchrotron optical emission are smaller than the 
shock radius, which indicates the synchrotron radio emis- 
sion and synchrotron optical emission are more diffuse. 
Actually, synchrotron optical emission from the cores of 
elliptical galaxies has been claimed to be observed by 
Hubble Space Telescope {HST) with a typical luminos- 
ity ^ 10^5 - 10^9 erg s-i Hz-i (e.g., iChiaberge etaH 
120001) . The detection of synchrotron X-ray emission from 
elliptical galaxy cores is on the edge of currently feasi- 
ble Chandra observations (e.g., iTananbaumI l200ll . and 
private communication) . 

One point which needs to be emphasized is that the 
mechanism for synchrotron radio emission discussed here 



is different from that produced by jets. Current under- 
standing of the AGN phenomena divides the activity 
into two phases: the 'optical mode', immediately after 
the central outburst, and 'radio mode' at low Eddington 
ratios between major outbursts, when powerful jets are 
emitted. In particular, here we have discussed the par- 
ticle acceleration during the optical mode, which should 
exist for each quasar event associated with the outburst 
phase. Though the synchrotron radio emission from large 
scale jets may be stronger than the presented mechanism, 
our mechanism is also important for the synchrotron 
emission at the central part of elliptical galaxies and it 
could also be important for the evolution of the galaxy 
since the energy is released within the gaseous envelope of 
the galaxy. Definite observation of non-thermal emission 
from AGN hosting galaxy cores would strengthen the 
case for AGN feedback. Though there may be unresolved 
jets in the galaxy cores as pointed out by some observers 
(e.g., Ho 2008 and references therein), our mechanism 
contributes comparable emission to these small jets, if 
not larger. 
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